Thermoelectrical manipulation of nano-magnets 
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We propose a device that can operate as a magneto-resistive switch or oscillator. The device 
is based on a spin-thermo-electronic control of the exchange coupling of two strong ferromagnets 
through a weakly ferromagnetic spacer. We show that the local Joule heating due to a high con- 
centration of current in a magnetic point contact or a nanopillar can be used to reversibly drive the 
weak ferromagnet through its Curie point and thereby exchange- decouple the strongly ferromag- 
netic layers, which have an antiparallel ground state. Such a spin-thermionic parallel-to-antiparallel 
switching causes magnetoresistance oscillations where the frequency can be controlled by proper 
biasing from essentially DC to GHz. 



Manipulation of the magnetic state on the nanome- 
ter scale is the central problem of applied magneto- 
electronics. The torque effect^, which is based on the 
exchange interaction of electrons injected into a ferro- 
magnetic region, is one of the key phenomena leading 
to current-induced magnetic switching. Current-induced 
precession and switching of the orientation of magnetic 
moments due to this effect has been observed in many 
experiments^^&^ii^ 3 .. 

The efficiency of current-induced switching is, however, 
limited by the necessity to work with high current densi- 
ties. A natural solution to this problem is to use electrical 
point contacts (PCs). This is because now the current 
density is high only near the PC, where it can reac h 14 ! 15 
values of order 10 9 A/cm 2 . The characteristic energies 
supplied to the electronic system are determined by the 
voltage drop V, and are for V ~ 0.1 -j- 1 V (which is eas- 
ily reached in experiments) comparable to the exchange 
energy in magnetic materials. 

One of the characteristic features of such high current- 
density states is the possibility to create local heating 
in regions with a size of a few nanometers. The energy 
supplied to the electronic system in this restricted region 
results in an enhanced local temperature, which to a high 
degree of accuracy can be controlled by the applied volt- 
age. Electrical manipulation of nanomagnetic conduc- 
tors by such controlled Joule heating is a new principle 
for current-induced magnetic switching. In this Letter 
we discuss one of the possibilities of the thermoelectri- 
cal magnetic switching caused by a non-linear interaction 
between the spin-dependent electrical transport and the 
magnetic sub-system of the conductor due to the Joule 
heating effect. We predict that a special design of mag- 
netic PC can provide both a fast switching and a smooth 
change of the magnetization direction in nano-size re- 
gions of the magnetic material controlled by the voltage. 
We also predict temporal oscillations of the magnetiza- 
tion direction in such regions (which are accompanied 
by electrical oscillations) under an applied DC voltage. 
These phenomena are potentially useful for microelec- 
tronic applications such as memory devices and voltage 



controlled oscillators. 

Equilibrium magnetization distribution in a magnetic 
stack. The system under consideration is shown in Fig [I] 
where three ferromagnetic layers are coupled to a non- 
magnetic conductor. We assume the following conditions 
to be satisfied: The Curie temperature Tc of region 1 
is smaller than the Curie temperatures Tc°'' 2 ^ of regions 

and 2; in region 2 there is a magnetic field directed 
opposite to the magnetization of the region, which can 
be an external field, the fringing field from layer 1, or a 
combination of the two. We require this combined mag- 
netostatic field to be weak enough so that at low temper- 
atures the magnetization of layer 2 is kept parallel to the 
magnetization of layer due to the exchange interaction 
between them via region 1 (we assume the magnetization 
direction of layer to be fixed) . In the absence of an ex- 
ternal field and if the temperature T > Tc , the proposed 
F/f(N)/F tri-layer is similar to the antiparallel spin-flop 
'free layers' widely used in memory applications^. 

As T approaches T c the magnetic moment of layer 

1 decreases and the exchange coupling between layers 
and 2 weakens. This results in an inhomogeneous distri- 
bution of the stack magnetization. Euler's equation for 
the magnetization distribution that minimizes the free 
ener gy o f the system can be written as follows (see, e.g., 
Ref.E3): 



dx 
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Here the x- and z-axes are directed along the stack and 
the magnetization direction in the region 0, respectively; 
9{x) is the angle between the magnetic moment M(x) at 
point x and the z-axis; M = |M|; a ~ I/aM 2 , where 
/ ~ kT c and a is the exchange energy and the lattice 
spacing, respectively, k is the Boltzmann constant; j3 is 
the dimensionless constant of the anisotropy energy. Be- 
low we assume the lengths Li ; 2 of regions 1 and 2 to be 
smaller than the domain wall lengths in these regions. 

Using Eq. and the boundary condition at the inter- 
face S between magnetic region 2 and the non-magnetic 
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Figure 1: Orientation of the magnetic moments in a stack of 
three ferromagnetic layers (0, 1, 2); the presence of a magnetic 
field directed opposite to the magnetization is indicated by an 
arrow outside the stack. 
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Figure 2: Magnetic moment orientations in the stack at tem- 
peratures Ti < T < T c (1) . 

region 3 d6/dx\s = (see, e.g., Ref. fl8h one finds the 
magnetization in region 1 to be inhomogeneous 9±(x) = 
02x/Li, < x < L\, while the magnetic moments in 
region 2 are approximately parallel with the accuracy 
a 1 M?(T)/a 2 M 2 (T) < 1 (the subscripts 1 and 2 refer to 
the region 1 and 2, respectively); the dependence of their 
tilt angle 9 2 on H and T is given by the equation 



9 2 = D(H,T) sin 9 2 ; D(H,T) 



Li(L 2 gM 2 (T)) 
4ai Mf(T) 



(2) 



While writing Eq. © we neglected the anisotropy energy 
in comparison with the magnetic energy in region 2. 

At low temperatures (D(T,H) < 1) Eq.© has the 
single root 6 2 = and hence a parallel orientation of all 
the magnetic moments in the stack is thermodynamically 
stable. At temperatures T >T\ (here T\ is the tempera- 
ture at which D(T, H) = 1) one has D > 1 and a second 
root 9 2 7^ appears in addition to the root 9 2 — 0. The 
parallel magnetization corresponding to 9 2 = is now 
unstable and the magnetization direction of region 2 tilts 
with an increase of T (see Fig© until at T > T c (1) the 
exchange coupling between layers and 2 vanishes and 
their magnetic moments become antiparallel. 

Current-voltage characteristics of the stack under 
Joule heating. If the stack is Joule heated by the cur- 
rent J its temperature T(V) is determined by the heat- 
balance condition 



JV = Q(T), J = V/R{6 2 ) 



(3) 



and Eq. © that determines the temperature dependence 
of 9 2 (T(V)). Here Q(T) is the heat flux from the stack, 
R{9 2 ) is the stack resistance. Eqs. © and © define the 
current-voltage characteristics (IVC) of the stack J = 
R{9 2 {V))V in a parametric form. 

As follows from the previous paragraph, in the whole 
intervals T(V) < Ti and T(V) > T c (1) the stack re- 
sistances are R(0) and R(ir), respectively, that is IVC 
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Figure 3: Numerical results for the I — V characteristics 
(IVC) of the magnetic stack. Here R(02) ~ R+ — R- cos #2, 
R-/R+ = 0.2, D = 0.2; J c = V c /R+. The branches - a' 
and b — b' of the IVC correspond to the parallel and antipar- 
allel orientations of the stack magnetization, respectively; the 
branch a~b corresponds to the inhomogeneous magnetization 
distribution shown in Fig[2] 



branches J = R(0)V and J = R(ir)V are linear at 

V < Vi = Vfl(0)Q(Ti) and V > V c = sj R{-k)Q{tP), 
respectively. At Vi < V < V c the stack temperature is 
Ti < T(V) < and the magnetization direction of 

region 2 changes with a change of V and hence IVC is 
non-linear there. Differentiating Eqs. (|2l3p with respect 
to V one finds that inside this interval the differential 
conductance Rd(V) — dJ/dV is negative if d[G{9 2 ){\ — 
D sin 9 2 /9 2 )}/d9 2 < (here D = D{H,T = 0)); e.g., for 
the stack resistance of the form R(9 2 ) = R + — i?_ cos 9 2 ) 
(here R± = (R(ir) ± i?(0))/2) one has dJ/dV < if 
D (H) < 3r/(l + 3r) where r = R-/R+, and hence IVC 
of the stack is N-shaped as shown in FigJ3] 

Self- exciting electric, thermal and magnetic direction 
oscillations. The set of equations describing thermal and 
electrical processes in the stack connected in series with 
an inductance C and DC voltage drop V (see FigUJ) is 



Cv^ = J 2 R(9 2 )-Q(T); 



£^L + JR(e 2 ) = V (4) 



where Cy is the heat capacity. For the case that the mag- 
netic moment relaxation is the fastest process the tem- 
perature dependence of 9 2 = 9 2 (T(t)) is given by Eq.©. 

The set of equation Eq.© has always a steady-state 
solution Eqs. (|2l3p . A study of its stability with respect to 
small perturbations shows that inside the interval Vi < 
V < V c this solution is unstable if C > C cr = t v \Rd(V)\ 
where t v ~ (TCyR+ /V 2 ) is the characteristic time of the 
voltage evolution and Rd = dVj d J < is the stack differ- 
ential resistance. As a result periodic oscillations of the 
current J(i) and the voltage drop on the stack V(t) spon- 
taneous arise: the current J(t) and the voltage drop on 
the stack V(t) periodically move along the limiting cycle 
in the plane (J, V) as is shown in FigfSl The temper- 
ature of the stack T(t) and the magnetization direction 
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Figure 4: Equivalent circuit for a Joule-heated magnetic 
stack. A differential resistance Rd{V) biased by a fixed DC 
voltage V is connected in series with an inductance £; V(t) is 
the voltage drop over the stack and J(t) is the total current 




Figure 5: Numerical results showing spontaneous oscillations 
of the current J(t) and the voltage drop V(t) over the stack for 
the parameters R-/R+ — 0.2, Do = 0.1 and [C — C cr )/jC cr = 
0.25; J c = V c /R+. The time development of the current and 
the voltage drop over the stack follows the dashed line or 
the thin solid line towards the limiting cycle (thick solid line) 
depending on the initial state. 



6*2 (i) adiabatically follows these electric oscillations ac- 
cording to the following equalities: Q(T(t)) — V(t)J(t) 
and 6 2 {t) = 9 2 {T{t)) (see Eq.©) For (C - C cr )/C cr <C 1 
the oscillation frequency is uj = 1 jt v . 

With a further increase of the inductance the character 
of the periodic motion changes its character: at C 3> 



C cr the current and the voltage slowly move along the 
branches of the IVC 0a and bb with the velocity J/ J « 
1/tj (here tj = C/R+ is the characteristic time of the 
current evolution) fast switching between them at points 
a and b with the velocity ~ l/t v . Therefore, in this case 
the stack periodically switches between the parallel and 
antiparallel magnetic states. 

In conclusion, we have shown that Joule heating of 
the magnetic stack presented in FigfT] allows electrical 
manipulations of the mutual magnetization orientation 
of layers and 2 as follows: 

1) In the regime of a controlled current there is a hys- 
teresis loop: with an increase of the current along the 
a-a' IVC branch the magnetic directions are parallel up 
to point a. At this point the stack is heated up to the 
temperature T\ at which D = 1 and the parallel ori- 
ented state looses its stability. As a result the system 
switches to the b - b' branch of IVC at which T > T c (1) 
and the magnetic moment M 2 switches to the antiparal- 
lel orientation. With a further change of the current the 
stack remains on the b-b' branch of IVC until the current 
reaches point b where the stack is cooled below Tc and 
system switches again to the a — a' branch at which the 
magnetic moments are parallel. 

2) At the regime of a controlled voltage and a low 
inductance one may smoothly re-orient the magnetiza- 
tion directions from the parallel to the antiparallel (see 
FigEl) changing the applied voltage inside the interval 
Vi < V < V c that is moving along a — b branch of IVC. 

3) An increase of the inductance of the circuit under a 
fixed DC voltage (see Fig|4| allows to excite electrical pe- 
riodic oscillations accompanied by periodical switchings 
of the magnetization directions in the stack from parallel 
to the antiparallel orientation. 
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